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SECTION 1 


INTRODUCTION 


The introduction of closed-loop electrohydraulic fatigue machines in 
the last few years has made it practical to perform fatigue testing without 
restrictions on the order of loading or the number of consecutive load cycles 
of one amplitude. In fact, complete flight-by-flight load simulations are 


possible, 


Schijve in Reference 1 presents a rather complete discussion on the 
subject of aircraft structure fatigue. The discussion covers the physical 
and empirical limitations associated with using a linear damage theory and 
supporting data for life assessment problems involving a variable amplitude 
load sequence. He noted that flight simulation testing was essential for 


accurate assessment of fatigue life and crack propagation rate, 


The present effort is concerned with establishing methodology for 
developing simulation sequences. The approach taken included the develop- 
ment of strain based sequences from strain gage data recorded during 
bomber and fighter aircraft flight operations. During the development of 
these sequences, the strain gage data were edited to remove small cycles 
(less than 3800 psi peak to peak) and inactive portions of the flight, but the 
order in which the remaining peaks occurred during a flight was preserved, 
The stress magnitudes were expressed in terms of stress relative to zero, 
in the same manner as the data were recorded, so that there was no 
requirement to identify the mean and alternating parts of the cycles for 


these sequences, 


The strain gage data was also reduced to obtain statistical data on 
parameters considered necessary to describe the intensity of the flight-by- 
flight data. Statistical distributions were obtained for preflight stress, 


ground and flight l-g stresses, stress cycles per ground operation, stress 
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cycles per flight, flight and ground peak stresses, Different counting 
techniques were used to obtain three sets of distributions of peak stresses 
resulting from gust, maneuver and taxi induced loadings, These data were 
used to develop simulated flight-by-flight sequences. It should be noted 
that G.M. vanDijk in Reference 2 presents several methods whereby one 
can process flight loads data to obtain statistical parameters that can be 


used to generate simulated flight data. 


The strain based and simulated sequences were applied to coupon 
specimens, and the resulting fatigue lives compa red. Techniques were 
selected as a result of these comparisons for simulating bomber and 
fighter aircraft operational load histories. Additional strain based and 
simulation sequences were developed from alternate bomber and fighter 
strain gage data as a means of demonstrating the general applicability of 
the selected simulation techniques for a range of operational flight load 


intensities. 


Volume I of this report provides a description of the data and pro- 
cedures used to develop the flight-by-flight sequences, the experimental 
program and a discussion of the results, Volume Ii provides a descrip- 


tion of the computer programs used in generating the simulation sequences, 
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SECTION 2 


RECORDED DATA 


2.1 B-58 DATA 


The University of Dayton during the period 1965 to 1969 had frequency 
modulated (FM) magnetic tape recorders installed on 4 of the USAF B-58 
bomber aircraft. These recorders were used to collect flight loads data 


from the aircraft as they were used in the normal operating manner, 


The recording program collected time history data on airspeed, alti- 
tude, center-of-gravity acceleration, and eight strain gage locations. 
Several of the strain gage locations were duplicated on the left and right side 
of the aircraft. For the present effort the strain gage located on the lower 
side of the spar at wing station 438.3 and span station 56.5 was chosen for 
the source of data. The stress at this location is sensitive to wing bending 
loading and is negative when the aircraft is on the ground and is positive when 
the aircraft is in flight. The time-history strain gage data was processed 
in terms of stress relative to zero stress using a calibrated value of pre- 
flight stress, which varies as a function of the aircraft gross weight, as the 
initial stress value when the aircraft is parked on the ground, All B-58 
flights were designated as operational training flights. The mission pro- 
files included the normal mission segments of taxi, take-off, ascent, cruise, 
descent, and landing as well as special mission segments of supersonic 
dash, refueling, and low level. The low level mission segments included 
flights over specified low level courses at constant pressure altitudes and 


at fixed altitudes above the terrain. 


It was decided, based on the mean stress level and the Ko of the test 
specimen, that the B-58 data could be truncated so that minimum range 
between one reversal and the next could be 3800 psi. The flight recordings 


were then processed such that the peak stress of any cycle and the minimum 
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stress of the cycle were written on digital tape. The position of a switch on 
the landing gear also was used to identify whether the stress peak occurred 
when the aircraft was on the ground or in-flight, As each flight was pro- 
cessed the statistical variance of the stress peaks and valleys was calculated 


by the equation: 


N /N 2 
N 5 7 - 5 
r® i 2 i 
| ; _ i=l i= 
Stress Variance = —ii #~#~€~€|CM (1) 


where S is the value of the stress for a peak or valley. 


One of the objectives of the project was to investigate the effect that 
a mild or severe load history might have on the accuracy of a simulation 
technique; therefore the flights were divided into two sets: a low intensity 
set and a high intensity set. The variance was used as the parameter to 
assign a particular flight to either the high or the low intensity set. A value 
for the variance of the airborne segment of a flight of 30.13 x 10° (wet)* was 
selected as the dividing point. Using this criteria the B-58 data set was 
divided so that 343 flights were assigned to the low intensity set, defined as 
Sequence A, and 282 flights were assigned to the high intensity set, defined 


as Sequence B, 


The low intensity data contained 197,448 stress reversals or 98,724 
complete cycles. Of these, 2097 were taxicycles. The total data set con- 
tained 194, 287 reversals which defined primary peaks (a peak between 
successive mean crossings) and a total of 3160 secondary reversals that 
were superimposed on the primary reversals. The average of the l-g trim 
stress during the airborne portion of the flight was 12,920 psi with a 
variance of 12.1x 10° feat”: This average is the average of the l-g 
stresses at the time of a stress reversal for the 193, 255 peaks, see Table l. 
The high intensity data set contained 86,121 stress reversals or 43,060 


complete cycles, and 3226 of these cycles were associated with taxi or 
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ground operations. A total of 2499 reversals were secondary reversals. 
Table 1 is a summary of the statistical properties of the low and high 


intensity data sets for the B-58, 


An examination of the table reveals some numbers that require addi- 
tional explanation. As stated early the dividing line between the low and 
high intensity flights was based on a value of the variance of 30.13 x 16 (ues, 
That number is the variance of the stress reversals about the mean for that 
particular flight. The variance listed in Table 1 is the variance of all 
primary peaks about the mean of the entire data set (i.e., 190,438 points). 
When the variance is calculated for the stress reversals that are primary 
peaks and the deviation is measured from the instantaneous l-g trim stress to 
the stress peak (or valley) then the variance for the Sequence A set is 
19.8 jews and for the Sequence B set is 24.9 aa”. These latter values 


labeled the 1,0G Variance were calculated as follows, 


MN; 
>, 607" 
:; — ha 
1.06 ees = 2 
‘ lance = 7G (2) 
) N.} - 1 
1 
i={ 
where 
.th | | : Ah i... 
5. = value of } peak or valley stress inthe i flight 
N. = number of stress peaks for flight segment (ground or airborne) 
for it flight 
M = number of flights in data set 
G. = l-g trim stress for flight i 


i 
Also some comments are required about the means, The mean of the 
primary peaks is just the simple average of the 190,438 peaks. The mean 
of the weighted 1l-g flight stress is the average of the l-g trim stress at the 
time of a stress reversal. The l-g trim stress values were determined by 


measuring the average l-g stress value just preceding each stress reversal 
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which yielded a time history of the 1-g trim stresses for the entire flight. 
The fact that the averages of the primary peaks and the weighted l-g stress 
are essentially the same is an indication that the deviations of the stress 
from the 1-g value is the same for increasing and decreasing deviations 
(i.e., plus and minus incremental stresses), This does not imply that the 
individual cycles are symmetrical about the l-g trim value but only that 
when the total data set is examined there is symmetry between the positive 
and negative peaks. The fact that the l-g stress on the ground is different 
from the preflight stress indicates that some of the data are from landings 
and some are from the take-off run just prior to lift off. The fact that the 
mean of the ground data and the mean of the 1-g ground data are different 


indicates that there are larger negative peaks than positive peaks, 


The preflight stress is the true magnitude of the wing stress when the 
aircraft is parked with engines off. These values were used as the starting 


point for processing of the strain gage data. 


Test sequences were also developed wherein taxi loads were omitted; 
however one ground load equal to the preflight stress was inserted at the 
beginning of each flight. The Sequence A data without taxi loads is identified 
as Sequence C, and the Sequence B data without taxi loads is identified as 


Sequence D, 


2.2 F-106 DATA 


The University of Dayton during the period from August 1970 to 
September 1971 had digital tape recorders installed on four F-106 Inter- 
ceptors. In addition to the usual VGH data, strain gage data was recorded 
from a gage on the lower shoulder of spar number four at wing station 472 
at buttock line 44. 25. The recordings began at engine startup and ended at 
engine shut down, and therefore include taxi, take-off, and landing as well 


as the airborne portion of the flight. 
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The strain gage location was inboard of the landing gear and therefore 
the stress was tension when the aircraft was on the ground as well as when 
it was inthe air. Consequently the strain at this location was not charac- 
terized by the typical ground-air-ground stress reversal. In fact the stress 


on the ground and the l-g stress in-flight were the same. 


The data was divided into two sets: a low intensity and a high intensity 
set. The variance, Equation (1), of the stress history was used as the 
measure of severity. The dividing point between low and high was chosen 


as 19.36 x 10° sii. 


The low intensity data set had 381 flights with 52,128 stress reversals; 
51,462 of these occurred during the airborne portion of the flight, and 
48,446 were primary peaks, The average number of primary peaks per 
flight was 127, The variance, Equation (2), of the airborne primary peaks 


| | Z 
about the l-g trim stress was 11.18 (ksi) . 


The high intensity data set had 236 flights with 39, 752 stress reversals; 
32, 636 of the reversals were airborne primary peaks. The average number 
of stress reversals was 138 and hence there was an average of 76 complete 
cycles per flight. The variance, Equation (2), of the peaks and valleys about 


; 2 
the l1-g trim stress was 28,2 (ksi) . 


Table 2 is a summary of some of the more significant facts about the 
data sets. Sequence E was the low intensity set and Sequence F was the high 


intensity set. 
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SECTION 3 


TEST SPECIMEN 


The test specimen was machined from 7075-T651 aluminum plate. 
The specimens were machined to the dimensions shown in Figure 1, The 
stress riser hole in the center of the specimen was first drilled to 0. 375- 
inch diameter on a numerically controlled vertical milling machine. It was 
later redrilled to 0. 386-inch diameter on a bench drill press. The speci- 
mens were machined from 1/4 inch rolled plate without any additional 
machining on the flat faces. In order to eliminate any biasing of the test 
results due to specimen location within the two plates that were used, the 
specimens were assigned consecutive numbers from random locations in 
the plates. The specimen locations are shown in Table 3, The length is 
the roll direction of the plate which is the 11 inch dimension of the specimen, 
It was the intent of this effort to use the same specimen that was used on a 
previous program, That program, the results of which are reported in 
References 3 and 4, used a specimen that had a 0.375-inch drilled hole in 
the center. The specimen for the previous program had a fatigue life of 
50,500 cycles when tested at a constant amplitude of +5 KSI to +25 KSI net 
section stress. However, the specimens from the current program with a 
0. 375-inch diameter drilled hole had a fatigue life in excess of 250, 000 


cycles for the same loading. 


In an attempt to make the fatigue strength of the specimens from the 
two programs similar, the 0.375-inch hole was reamed to 0. 386-inch dia- 
meter, This reaming reduced the fatigue life to 67, 200 cycles for the con- 
stant amplitude loading. Since this life was still greater than the 50,500 
from the previous program, the following additional modifications were 
attempted: (1) the drilled 0. 375-inch diameter hole was stress relieved at 
200°F for 24 hours, yielding a fatigue life of 105, 200 cycles; (2) the 0. 375- 


inch hole was drilled out to 0.386-inch diameter, yielding a fatigue life of 
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58,200 cycles; and (3) the 0.375-inch hole was drilled out to 0.386-inch 
diameter and then stress relieved at 200°F for 24 hours, yielding a fatigue 
life of 61,100 cycles. Since the stress relieving after redrilling did not 
cause any additional reduction in fatigue life the specimens were just 
redrilled to 0.386-inch diameter without stress relieving. The final con- 
figuration had a fatigue life of 58,200 cycles when tested at a constant am pli- 


tude of +5 KSIto +25 KSI. 


The static material properties of the plates were determined and are 
presented in Table 4. Tensile specimens were machined from random loca- 
tions throughout the plates. The results indicate that the plates were uni- 


form in tensile properties. 
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SECTION 4 


TESTING PROCEDURE 


All cyclic testing was performed using the University's MTS closed- 
loop electrohydraulic fatigue test system which was programmed by signals 


from an FM recorder. 


4.1 GENERATION OF TEST SEQUENCE TAPES 


The loading spectra which resulted from the editing of the actual 
flight recordings was a digital tape of peak and valley stresses in the same 
sequence as they occurred during the flights. Now since 343 flights of 
Sequences A or C would not contain enough cycles to cause a failure, the 
343 flights were arranged in a second and third order to give a test sequence 
of 1029 flights. The order of the cycles within a given flight was not 
rearranged. After the 1029 flights were applied to the specimens the same 
sequence was started again. For the B and D sequences, the 282 flights 
were rearranged into 5 additional orders to give a test sequence of 1410 
flights. Again, the order of the cycles within a given flight was not changed. 
After the 1410 flights were applied to the specimens the sequence was 


restarted. 


The same procedure was used for the F-106 data. The 381 flights of 
Sequence E were rearranged 11 times to give a test sequence of 4572 flights. 
The 236 flights of the Sequence F data were rearranged 14 times to give a 
test sequence of 3540 flights. As with the B-58 data the order of the cycles 


within a given flight was not changed during the rearrangement process. 


As each flight of the B-58 data set was copied from the original edit 
files, the preflight stress and two calibration signals were inserted into 


the data. The preflight stress was inserted at the beginning of the ground 
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data and the calibrations signals between the ground and airborne data. For 
spectrums C and D which did not contain taxi cycles the preflight stress 
was followed immediately by the two calibration signals. Since the F-106 
data did not have a negative strain for ground operations the calibration 
signals were inserted into the data differently. The preflight stress of 
+8360 psi was written 4 times at the beginning of each flight, a signal equal 
to +14, 960 psi was written 4 times in succession once per flight if and when 
the stress peak following the calibration value exceeded 16,060 psi. The 
second calibration signal was inserted this way so that the calibration signal 
would not introduce an extra stress reversal. If a flight did not have a peak 
greater than 16,060 psi then the second calibration voltage was omitted for 


the flight. 


After the digital tape was created it was used as an input to the Wright- 
Patterson Air Force Base Analog/Hybrid computer. The analog signal was 
generated by fitting half sine waves between the peaks and valleys. The 
output of the analog/hybrid computer which was in an analog voltage form 
was recorded on a frequency modulated (FM) tape recorder. The sinusoidal 
frequency during the digital to analog conversion was 80 Hz, however, when 
the FM tape was played back into the MTS testing system the tape was slowed 


down so that the testing frequency was 5 Hz. 


Oscillograph playbacks of typical flights are presented in Figure 2. 
Figure 2(a) and 2(b) are the same B-58 flights with and without taxi cycles. 


Figure 2(c) is a section from an F-106 spectrum. 


In an attempt to control the amount of plasticity at the edge of the 
stress riser hole the net section stress was limited to 30, 000 psi tension. 
With a stress concentration factor of 2.54 on the net section stress this 
limited the maximum stress to 76,200 psi. The 0.2% offset yield stress was 


75,600 psi. The net section compression stress was limited to -15,000 psi. 
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4.2 SPECIMEN INSTALLATION 


The specimens were installed in the testing machines using the gripping 
device shown in Figure 3. During installation the bolts were tightened first 
while the specimen was loaded slightly in tension (100 pounds) then the wedges 
were tightened. During testing the tensile loads were transfered to the speci- 
men by shear loads on the bolts and the compressive loads were transfered 
directly to the specimen by the wedges. The use of the wedges in combination 
with the bolts allowed tension-compression testing without any backlash in 
the bolt holes. Using the relatively stiff fixture and specimen it was possible 


to test the specimens without using an anti-buckling fixture. 


4.3 CYCLIC TESTING 


Cyclic testing was accomplished by using the FM tape signal as the 
command signal for the closed-loop testing system. The feedback signal for 


control was from a strain gage type load cell. 


The FM recorder/reproducer used was a 14 channel system with two 
channels used for WOW and flutter compensation so 12 channels were used 
for test sequence data. Each test sequence (A, B, C, D, E, and F) was 
limited to approximately 300,000 cycles before the sequence was repeated. 
At a tape speed of 1.875 inches/second it required three channels to record 
the 300,000 cycles on 10.5 inch diameter reels. Each channel (track) of the 
tape contained approximately 100,000 cycles and provided 6 hours of testing 
time. Sequence A for example started on track 1] continued on track 2 and 
was completed ontrack 3. Sequence C started on track 4 continued on 
track 5 and was completed on track 6. Each track had a section of zero load 
recorded at the beginning and end of the test sequence. During testing the 
reproduce recorder would be started, and while the signal was at the zero 
load level, the output would be connected to the MTS system. After a track 


was reproduced, the output would be disconnected from the MTS system 
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while the signal was at zero load, and the tape rewound. The next track 
would then be used to program the machines. If after tracks 1, 2 and 3, 
for example, had been used and the specimen had not failed, then testing 


was continued by using track 1] over again and then track 2 again and so on. 


Since the system used is completely analog, it is necessary to monitor 
the load/cell output to ensure that the gain of the FM system has not changed 
or that a DC offset has not occurred. The accuracy of the loadings was main- 
tained by monitoring on an oscilloscope the calibration signals which were 
recorded periodically on the tape. In order to obtain a high sensitivity on 
the oscilloscope, the load cell output was input into a chopped-beam oscillo- 
scope through two "Tektronic'' Model 5A13N Differential Comparator Ampli- 
fiers. In using these Differential Comparator Amplifiers, the lower cali- 
bration voltage was set on one amplifier and the higher calibration voltage 
was set on the other amplifier. When the test spectrum was being applied 
to the specimen and when the periodic calibration pulses occurred, the 
difference between the actual load cell voltage feedback and the required 
calibration voltages was displayed on the oscilloscope screen. If both the 
gain and the DC level were correct then the output during the two calibra- 
tion pulses would be zero. Since the output should be zero, the sensitivity 
(vertical deflection) of the oscilloscope could be set very high. The sen- 
sitivity was set so that 0.1] inches on the screen was equivalent to 30 psi 
net section stress. By frequent monitoring of the calibration pulses in the 
load cell output signal, the stress levels for the calibration signals were 
maintained within +60 psi. 

Cyclic testing was conducted on a 24 hour a day 5 days a week basis. 


The specimens were maintained at zero load over the weekends. 
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SECTION 5 


SIMULATED TEST SEQUENCES 


The strain gage data was processed to obtain statistical information 
on parameters considered appropriate to characterize the intensity of the 
flight data. Initially these processing efforts were restricted to the strain 
data that was used in the development of Sequences A and E as explained in 
Section 4. The strain data was processed to obtain statistical distributions 
of preflight stress, ground and airborne l-g stress, stress cycles per 
ground and airborne operations and peak stresses. It was contemplated 
that the information on stress peaks could significantly influence test life. 
Consequently, different counting methods were used to obtain alternate sets 
of peak stress distributions. Having obtained these data sequences were 
developed and applied to test specimens to simulate Sequences A and E, 
The stress peak counting method employed to develop the simulation 
sequences that produced test lives that correlated with the strain based 
Sequences A and E test lives was selected for additional evaluation. 
Sequences were then developed to simulate the strain based Sequences Bb, 


C, D, and. F. 


5.1 B-58 SIMULATIONS 
5.1.1 Folded Distribution 


The first method used to generate a simulated test spectrum 
was the folded distribution of combined positive and negative primary peaks. 
A primary peak is defined as the maximum (or minimum) value of stress 
occurring between two successive crossings of the l-g stress level by a 
continuous stress time-history trace. In order to generate the input data 
for the simulation program, the actual strain gage data was processed by a 
computer program to generate statistical distributions, in the form ofa 


cumulative frequency distribution, of incremental stress peaks and valleys. 
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In order to conserve testing time only those peaks (valleys) greater than 
+ 3000 psi were used, Thus, an incremental peak stress represents an 


alternating stress superimposed ona l|-g trim stress, 


Separate distributions were generated for the ground opera- 
tions and the airborne flight segments. The frequency distributions of 
incremental stress peaks for the B-58 Sequence A data set is presented in 


Table 5 for the ground operations and also for the airborne operations, 


Table 6 is a listing of the folded distribution (i,e., the sum of 
the negative and positive peaks) for the ground data for the B-58 Sequence A 
data. Since the parameter that is needed to generate a simulated test 
sequence is the relative probability, Table 6 lists the sum of the positive 
and negative peaks for the frequency and cumulative frequency distributions 
from Table 5 rather than the average of the number of peaks and valleys. 
The probability listed is the probability of a peak being equal to or greater 
than the absolute value of the stress value listed. Table 7 is a listing of 


the same type of data for the airborne flight segments. 


Folded distributions of the primary stress peaks were used to 
generate random test sequences, A flight within the sequence was assembled 


by the following sequential steps: 
(1) randomly select a preflight stress from the cumulative 
probability of preflight stress data (Table 8); 


(2) randomly select a ground 1l-g trim stress (Table 9); 


(3) randomly select the number of stress cycles for the 


ground operations for that flight (Table 10); 


(4) randomly select the magnitude of the cyclic stress for the 
first ground stress cycle (Table 6), add and subtract this value from the 
ground l-g trim stress selected in step (2), and continue to repeat step (4) 


until the number of ground cycles required by step (3) is obtained; 


ad 


(5) randomly select a l-g trim stress value from Table 11 


for the airborne segment of the flight; 


(6) randomly select the number of stress cycles for the 


airborne segment of that flight (Table 12); 


(7) randomly select from Table 7 the magnitude of the cyclic 
stress for the first stress peak of the airborne segment of the flight, add 
and subtract this value from the l-g trim stress selected in step (5) to form 
the first stress cycle, and repeat step (7) until the number of cycles selected 


in step (6) is obtained. 


A random number between 0. 00000000000 and 1, 00000000000 
was selected by a random number generator then used to interpolate in the 
appropriate probability table to select values for the parameters in steps (1) 


thru (7). 


This procedure completes the first flight. Now go back to 
step (1) and repeat the entire procedure until enough flights and cycles are 
generated to fill one track of the FMtape. The procedure is repeated again 


for a second and third FM channel. 


In the actual simulation of the B-58 Sequence A data the first 
channel contained 420 flights with 105,651 stress cycles, the second channel 
439 flights with 105, 393 stress cycles and the third channel 459 flights with 
105, 303 stress cycles. 


Table 13 is a listing of a typical summary sheet for the first 
50 flights. The summary sheet is self explanatory except for two items, 
i.e., the cumulative number of peaks and points. The number of peaks is 
the sum of the number of ground and airborne stress reversals (2 per cycle), 
the number of points is equal to the number of peaks plus the number of 
points used for the preflight stress (4) and the ground calibration pulse (4) 


and the airborne calibration pulse (4) so that each flight has 12 more points 
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than peaks except the first flight which has 8 points of zero stress at the 
beginning. A summary of the total data set on each of the three channels 
is presented in Table 14. Table 14 also contains the summary data for the 


entire data set. 


Note that there are two variances listed in Table 14 for both 
the ground and the airborne data. The one labeled stress variance was 
based on Equation (1) and the one labeled 1. 0G Variance was based on 


Equation (2), 


Table 15 is a listing of stress reversals for the start of each 
channel, This print out is from the digital tape which was used as input to 
the hybrid computer, A continuous analog signal for each channel was 
obtained by forcing a half cycle of a sine wave between the digital peaks 
and valleys. The resulting signals were then recorded on FM tape and later 


used to program the fatigue machines. 
5.1.2 Bivariate Distribution of Positive and Negative Primar Peaks 


The second method of generating a simulated test spectrum 
requires the strain gage data to be processed by a different computer pro- 
gram which sorts the primary peak data into a bivariate distribution of the 
peak stress in a cycle and the immediately following minimum stress. 

Table 16 is a listing of the maximum stress and the corresponding minimum 
stress of each cycle of the ground data for the B-58 Spectrum A data set. 

In this data set only those cycles wherein the maximum stress was more 
than 3000 psi greater than the l-g trim stress were retained. The 
immediately following minimum did not have to exceed a minimum threshold 
but could have any value provided it was less than the l-g trim stress. 
Table 17 is a bivariate distribution of peaks and corre sponding valleys for 


each cycle of airborne data, 


Table 18 is a listing of the cumulative probabilities for the 


peak stresses of the ground operations data and Table 19 is a listing of the 
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cumulative probabilities for the immediately following minimum stress for 


each peak stress for the ground operations data. 


Tables 20 and 21 are the same type of data for the airborne 
segment of the flights. 


All other data is the same as that used for the folded distribu- 


tion described in Paragraph 5.1.1. 


The method for generation of a test sequence from this data 
is the same as described in Paragraph 5.1.1 except for steps 4 and 7. For 
the bivariate distribution of peaks and valleys method, one must first select 
a random number and then find the corresponding magnitude of the peak 
stress (Table 18 for step 4 and Table 20 for step 7). This value of cyclic 
stress is added to the l-g trim stress obtained in steps 2 or 5. Next, with 
a second random number, one selects a minimum stress from the distribu- 
tion of valleys for that magnitude of peak stress (Table 19 or Table 21). 
The magnitude of the valley is then subtracted from the l-g trim stress to 


determine the extremes of a complete stress cycle. 
A summary of the first 50 flights is presented in Table Ze. 


A summary of the three channels of data used for the simula- 


tion of the B-58 data Sequence A is presented in Table a3. 


Table 24 is a listing of the start of the data as it appeared on 
the digital tape that was used as input for the Hybrid computer, 


5.1.3 Limited Folded Distribution 


The third method used to generate a simulated test spectrum 
was identified as the Limited Folded Distribution. This method is similar 
to the Folded Distribution (Paragraph 5.1.1) in that the cyclic stresses are 
selected from a folded distribution of the primary peaks and are applied 
symmetrically about the l-g trim stress; however, the other parameters 


(preflight stress, number of ground stress peaks per flight, ground l-g 
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stress, number of airborne stress peaks per flight and the l-g stress for 
the airborne flight segment) were selected from a small discrete distribu- 
tion, The input data had only a few values for these parameters and the 
program did not interpolate for intermediate values, Tables 25 through 31 
list the input data for the simulation program. Table 32 is a summary 
listing of the first 50 flights generated by the program, Note that the 
number of stress reversals for the ground segments have discrete values 
of 0, 4, 8, 12, and 78 and the airborne segments 74, 238, 454, 784, and 
1620, respectively. Also note that the preflight and l-g trim stress have 


only a few discrete values, 


5.2 F-106 STRESS SIMULATION 


5.2.1 Bivariate Distribution of Mean and Alternating Stresses 


The basic principle of this method is to form a stress cycle 
from mean and alternating stress values chosen randomly from a bivariate 
distribution of mean and alternating stresses. First, a mean stress value 
is randomly chosen from the cumulative distribution of mean stresses, 
Then an alternating stress value is randomly chosen from the cumulative 
distribution of alternating stresses for the chosen mean stress value. A 
stress cycle is then formed by adding the alternating stress to the mean 
stress thereby obtaining the peak value of the stress cycle and by sub- 
tracting the alternating stress value from the mean the minimum value of 
the stress cycle is obtained. Separate bivariate tables are used for the 
ground and airborne data. To complete the desc ription of a flight, the 
number of ground cycles and the number of airborne cycles for the flight 
are randomly chosen from the distributions of the number of ground and 
airborne cycles per flight. The input data for the computer program which 
generated the simulated flights for the F-106 Sequence E (Low Intensity) 
data set are presented in Tables 33 through 38. The cumulative distribu- 


tions are entered with random numbers which are equated to the probabilities, 
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Summary information for the first 50 flights is presented in 
Table 39. The variance, Equation (1), of the stress reversals about the 
mean of the reversals is less than the variance, Equation (2), of the 
reversals about the l-g trim value. The l-g trim stress was taken as a 
constant value of 8360 psi. The difference in the variance is due to the 
unsymmetrical nature of the loadings with the positive incremental stresses 
being larger than the negative incremental stresses, The mean of all of 


the peaks and valleys is 10,607 psi whereas the 1-g trim stress is 8360. 


Using the bivariate distribution of mean and alternating stress 
10, 094 flights were generated which contained a total of 570,528 stress 


reversals or half that number of cycles. 


The composite summary of all the flights is presented in 


Table 40, 


5.2.2 Separate Positive and Negative Distributions 


The second method used to generate simulated flights of the 
F-106 Sequence E data was identified as the Separate Positive and Negative 
Distributions. This method requires the random selection of a positive 
incremental stress and then a negative incremental stress from separate 
distributions of the primary peaks. The positive incremental stress is 
added to the l-g trim stress (8360 psi) to obtain a stress peak and the 
negative incremental stress is added to the 1l-g trim stress to obtain a 
stress minimum. Separate distributions are used to describe the ground 
and airborne data. To complete the description of a flight, the number of 
stress cycles for ground operations and the number of stress cycles for the 
airborne segment of the flight are selected randomly from distributions of 
the number of stress cycles for ground and airborne flight segments per 


flight. 


The cumulative probabilities for the incremental stresses for 


the F-106 Sequence E data are presented in Tables 41 through 44. The 
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cumulative distribution of the number of stress reversals per flight (not 
cycles) is presented in Table 45 for the ground operations and in Table 46 
for the airborne operations. A summary of the first 50 flights generated 
by this procedure is presented in Table 47. A summary of all the flights 
generated is presented in Table 48. There were 12,111 flights with 
560,301 stress reversals used to fill three channels of the FM tape. 


5.2.3 Limited Separate Positive and Negative Distributions 


This method is similar to the method presented in Paragraph 
5.2.2, The stress cycles are formed by randomly selecting a positive and 
a negative incremental stress from separate distributions and adding them 
to the constant 1-g trim stress 8360 psi. However, the number of cycles 
in the ground and airborne segments are chosen from a small discrete 


distribution, 


The input data for the computer program is presented in 
Tables 49 through 54. A summary of the first 50 flights is presented in 
Table 55, 


Table 56 is a presentation of the summary of all the flights 
for this method of simulation. A total of 12,495 flights with 552,579 stress 


reversals were generated, 


5.3 SIMULATIONS FOR DATA SETS B-58 SEQUENCE B, C, D AND 
F-106 SEQUENCE F 
After tests using the simulated sequences which were discussed 
in Paragraphs 5,1 and 5,2 were completed, the remaining data sequences 
were simulated using the best simulation technique for each type aircraft. 
For the B-58 that was the Bivariate Distribution of Positive and Negative 
Primary Peaks, For the F-106 the Separate Positive and Negative 


Distribution was used for the stress data. 
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B-58 Simulations for Sequence B, C, and D 


The B-58 Sequence B data which were the high intensity flights 
were used to generate a test sequence using the bivariate distribution of 
peaks and corresponding valleys. This procedure is described in Para- 
graph 5.1.2. The input data is presented in Tables 57 through 66. The 
summary data of the simulated flights is presented in Table 67 for the first 


50 flights and in Table 68 for the entire test sequence, 


The simulated test spectra for the Sequence C and D data, 
which is really the B-58 low and high intensity data sets but without taxi 
cyclic stresses, used the same input data as the A and B sequences. The 
results of the simulation for the first 50 flights for Sequence C simulation 
is presented in Table 69 and the summary for all the flights in Table 70. 
The results of the simulation for the first 50 flights for Sequence D simu- 
lation is presented in Table 71 and the summary for all the flights in 


Table 72. 


5,3,2 F-106 Simulations for Sequence F 


The F-106 high intensity data, Sequence F, was used to 
generate a sequence of test loads using the method presented in Paragraph 


5, Za Zs 


The input data which was obtained from processing the strain 


gage data into statistica] parameters is presented in Table 73 through 78. 


The summary of the first 50 flights is presented in Table 79 
and the summary data for all of the Sequence F flights in Table 80. 
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SECTION 6 


TEST RESULTS 


All of the testing was conducted in the University's structural test 
laboratory at a temperature of approximately 72°F. The test specimens 
were assigned to the various tests in a chronological and sequential order 


starting with specimen number 1 (see Table 3) for test number 1. 


In most instances six specimens were tested for each loading 
sequence. The final fracture was selected as the event to identify as the 
fatigue life of the specimen. The ratio of crack growth period to the total 
flights to failure varied from a minimum of approximately 0.07 (speci- 
mens with cracks growing on both sides of the starter hole) to a maximum 
of 0, 22 (specimens with a single crack on one side of the starter hole). 
The basic test data indicates that on the average 85% of the fatigue life was 


used to initiate a visible fatigue crack at the starter hole. 


6,1 FATIGUE LIFE 


The fatigue life of each individual specimen tested with the original 
B-58 data sequences is presented in Table 81. Table 82 is the results of 


the F-106 original sequences. 


The results of the fatigue tests using simulated test sequences are 
presented in Table 83 for the B-58 Sequence A and Table 84 for the F-106 
Sequence E. The comparison of fatigue lives must be made on the basis of 
flights to failure, not cycles, since the number of cycles per flight varies 


with the simulation method and the threshold level. 


Since the method which used the bivariate distribution of primary 
peaks and corresponding valleys (Paragraph 5.1.2) gave acceptable corre- 
lations for the B-58 Sequence A data, the limited distribution (Paragraph 


5.1.3) wa s not tested. The remaining sequences of the B-58 data set were 
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simulated using the method of bivariate distribution of peaks and corre- 
sponding valleys. The results of the tests using the simulated spectra for 


Sequences B, C, and D for the B-58 are presented in Table 85. 


The best correlation for the F-106 Sequence E data set was obtained 
for the limited separate peak and valley distributions; however, the regular 
separate peak and valley also gave acceptable correlation and since it was 
more complex (and therefore have broader application) we decided to use 
the regular separate peak and valley (Paragraph 5. 2, 2) distributions for the 
Sequence F simulations. The results for the F-106 simulations for 


Sequence F are presented in Table 86. 


The average life of the tests for each sequence is presented in 
Table 87. Only the flights to failure are listed in the table since the 
number of cycles per flight varied depending on the simulation method and 
the truncation level, The comparison must be made only on the flights 


since this was the only criterion used in the simulations. 


6,2 CRACK GROWTH 


The length of the fatigue crack was measured on most of the speci- 
mens. The measurements were made with a 30x telescope mounted ona 


micrometer slide as shown in Figure 4. 


The data presented in Figures 5 through 19 show the relationship of 
crack length versus flights for different test sequences after an initial crack 
length of 0.3 inch was established. A baseline crack length of 0.3 inch 
(distance from center of starter hole to crack tip) was selected to provide 
a common starting point since the start of visible crack growth occurred 
at different numbers of flights and cycles for the various test sequences, 
Crack length data is plotted versus flights rather than cycles since the 


criteria for the threshold stresses were different for each simulation 
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technique and for the original data. The plots were terminated when one 
side of the specimen fractured or when a visible crack was established on 


the opposite side of the starter hole. 
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SECTION 7 


DISCUSSION 


Perhaps the best way to visualize the objective and results of this 
program is to think of the program in two parts, the B-58 data and the 
F-106 data. 


7.1 B-58 DATA CORRELATION 


The results of all the testing using the B-58 data is presented pic- 
torially in Figure 20. These results, also see Table 81, indicate that the 
average number of flights to failure was approximately the same when 
using the Sequence A or Sequence B spectra. However, Tables 20 and 
65, which present the distribution of airborne peaks for Sequence A and 
Sequence B respectively, show that Sequence A contained over twice the 
number of peaks in comparison with Sequence B. These tables also show 
that approximately 50% of the peaks were less than 4,500 psi and that both 
data sets contained the same number of peaks above 8,000 psi. It would 
appear therefore that there was little difference in the severity of the 
Sequence A and Sequence B data sets and that the damage per cycle below 
stress levels of 5,000 psi was small or non-existant, It would also appear 
that the use of threshold values below 4,000 psi, in an effort to reduce the 


testing time requirements, did not significantly affect the fatigue lives. 


It should be noted that the separation of the original data into low 
intensity (Sequence A) and high intensity (Sequence B) data sets was accom- 
plished to demonstrate that the subsequent data simulation methods would 
be applicable when using different data sets. However, the separation of 
the B-58 data which is a more gust sensitive aircraft yielded data sets of 
approximately the same intensity (on a damage per flight basis) except that 


the low intensity data set contained a significantly larger number of low 
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value stress peaks. One could also conclude that the operational missions 
of the B-58 fleet were relatively consistent and that differences in resulting 


spectra were due to environmental atmospheric conditions. 


Referring again to Figure 20 one can see that when the taxi cycles 
were omitted from the Sequence A original data that the fatigue life in- 
creased from 2, 536 flights to 3,697 flights. This result was somewhat 
surprising at first; however, when one considers that when the taxi cycles 
were omitted the minimum stress on the ground was then the preflight 
stress and not as negative as the taxi cycles then it is not as surprising. 
Figure 2(a) and 2(b) depict the difference in the minimum ground stress 
very dramatically when the taxi cycles are omitted. We now believe that 
rather than using the preflight stress as the ground stress, if one is 
omitting taxi cycles, we should use the minimum stress in the ground 
mission segment (for that flight) as the ground stress. Reference 5 pre- 
sents some data on this for crack growth but not for crack initiation. For 
simulations we suggest that one use a distribution of minimum stresses 
for the ground data rather than a distribution of preflight stresses. In 
other words, the stress that was used as the ground stress for the 
Sequence C simulation which is listed as preflight stress (Column 2) in 
Table 29 should be replaced by the minimum ground stress data (Column 7) 
of Table 22. Additional insight to the effect of ground cyclic stresses will 


be obtained as a result of a current testing program (AF 33615-74-C-3007). 


The simulated D sequence does not follow the other trends, i.e., for 
all but the simulated D sequence the simulated sequences using the bivariate 
distribution of peaks and corresponding valleys have slightly shorter lives 
than the original strain gage data. No reason for this opposite result is 


known other than possibly experimental scatter. 


The tests using the folded primary peak spectra yielded a much 
shorter fatigue life and as expected is a more damaging sequence than the 


real flight situation, 
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A general comment is in order about the maximum stress in the 
sequence, As stated before, the maximum net section stress was limited 
to 30,000 psi; however, at the start of testing we used a truncation level 
of 34,500 psi and also 33,000 psi. The Sequence A spectrum was used for 
testing and resulted in longer lives for some of the specimens but also much 
larger scatter in the fatigue lives. The use of a spectrum that has tensile 
stresses that induce notch yielding appears to cause residual compressive 
stresses that cause longer fatigue lives. This phenomenon of retarded 
crack growth rate and increased fatigue life due to tensile overloads has 
been well documented in the literature (see, for example, References 5, 

6, 7, and 8). When one is using a test sequence that has stresses which 
could cause yielding at the notch it is very important not to exceed the 
maximum stress, even the normal scatter in the yield stress of the material 


may cause yielding in some specimens and not in others. 


We believe that the results presented in Figure 20 for the B-58 are 
a correct representation of the fatigue life of the specimens. The best and 
also good correlation was obtained using a simulated test sequence wherein 
the loads were randomly selected from a bivariate distribution of positive 


primary peaks and corresponding primary valleys. 


7.2 F-106 DATA 


All of the F-106 data set test results are presented pictorially in 
Figure 21. The results clearly show the effect of the severity of the two 
sequences, The E sequence being more mild than the F, Table 44 which 
is the distribution of primary stress peaks for the E sequence lists 962 
stress peaks greater than 10,000 psi whereas Table 76 lists 2419 peaks 
greater than 10,000 psifor the F sequence, Looking at the cumulative 
probability one can see that in the E sequence only 6.8% of the primary 
peaks are greater than 10,000 psi whereas for the F sequence 22,4% are 


greater than 10,000 psi. 
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The use of the bivariate distribution of mean and alternating stresses 
cause a much shorter fatigue life. This was due mostly to the shift in the 
mean stress from cycle to cycle during the simulation, This result is 
illustrated in Figure 22. The variance, Equation (2), of the incremental 
peaks and valleys of the cycles for the simulation using bivariate distri- 
butions of mean and alternating stresses was 28,5 KSI“ as compared to 
the variance of 22,6 KsI* for the simulation using the separate positive and 
negative distributions of primary peaks, Both variances were calculated 


using a constant l-g trim stress for 8360 psi. 


The simulation tests using separate distributions of primary peaks 
and valleys resulted in fatigue lives that were only slightly greater than 
those using the original strain gage data. The difference in the average 
life is about 11%; however, there is overlap in the scatter bands as can be 
seen on Figure 21, The fact that the separate peak and valley and the 
limited separate peak and valley gave the same fatigue life is not surprising 


since the F-106 data did not have a variable 1-g trim stress or a GAG cycle. 


7.3 CRACK GROWTH CURVES 


The crack growth curves presented in Figures 5 through 19 are 
evidence that the rate of crack growth is dependent on the loading spectrum. 
The data also indicates that after the crack has propagated to the point 
where the plots were started there are very few areas where crack growth 


was arrested for a time. 
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SECTION 8 


CONCLUSIONS 


The use of a load sequence wherein the stress reversals were 
arranged by a random selection from a bivariate distribution of the 
primary peaks and the corresponding valleys yielded fatigue lives that 
were equivalent to the lives using the original strain gage data for a 


gust sensitive aircraft. 


The use of a load sequence wherein the stress reversals were 
arranged by a random selection from separate distributions of the 
primary peaks and valleys yielded fatigue lives that were equivalent 
to the lives using the original strain gage data for a maneuver sensi- 


tive aircratit. 


The proper representation of the minimum stress in a flight is 
necessary for good simulation. The minimum stress for a good 
simulation, in which ground cycles are omitted, is the minimum value 


of the ground stress and not the pre-flight trim stress. 


Test sequences can be developed that provide proper simulation of 
strain history data. In addition to the distribution called for in 1 and 2 
above, statistical distributions of the number of stress cycles per 
flight for ground and airborne operations and a distribution of l-g 


trim stresses are required. 
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COVER PLATE 
SPECIMEN 


BRASS SHIM 


DRAW BOLT 


WEDGE 


Figure 3. Specimen Installation Procedure. 
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Telescope Setup for Crack-Growth Measurements. 


Figure 4. 
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Test Results for All B-58 Data Sets. 
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SUMMARY OF LOAD SEQUENCE FOR B-58 DATA 


ORIGINAL STRAIN GAGE DATA 


PRIMARY PEAKS 
Airborne Data 
No. of Peaks 
Mean (psi) 
Variance (psi)* 


Ground Data 
No. of Peaks 
Mean (psi) 
Variance (psi) 


WEIGHTED 1l-g FLIGHT STRESS 


Airborne Data 
No. of Peaks 
Mean (psi) 
Variance (psi) 

Ground Data 
No. of Peaks 
Mean (psi) 


PREFLIGHT STRESS 
Mean (psi) 
Variance (psi) 


Mean of the No. of Airborne 
Primary Peaks Per Flight 


Variance of Above 


Mean of the No. of Ground 
Primary Peaks Per Flight 


Variance of Above 


SEQUENCE A 
Low Intensity 
343 flights 


190, 438 
12,905 
33,52 =e 10 


3849 
-8047 6 
43.2x 10 


193,255 
12,920 
i2.1x 10° 


4193 
-6640 


-8647 
4. 3x10* 


655,22 
43.6x107 


1L..2) 
447.5 


SEQUENCE B 
High Intensity 
282 flights 


80,727 
14,780 


41.8 x 10° 


2896 
~8232 ; 
46.9 x 10° 


82,895 
14,690 


17'S Ox 10° 


3226 
-6373 


-8729 
2.5x10% 


286.27 


11.2x16 


LO, 31 
69.7 


TABLE 2 


SUMMARY OF LOAD SEQUENCE FOR F-106 DATA 
ORIGINAL STRAIN GAGE DATA 


SEQUENCE E SEQUENCE F 

Low Intensity High Intensity 
Number of Flights 38] 236 
Total No. of Stress Reversals 52,128 39,752 
Airborne Primary Peaks 48,446 32,636 
Variance of Primary Peaks LA. 18(ksi)@ 28.20(ksi)* 
Stress Reversals (Ground Operations) 666 652 
l-g Trim Stress (Constant) 8360psi 8360psi 

ky, iy 
2 
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TABLE 3 
SPECIMEN ASSIGNMENT CHART 


Specimen Specimen Specimen 


Specimen Location Specimen Location Specimen Location 
No. * Code ** No. * Code ** No, * Code ** 
] Eega | 34 D 4-4 | 67 CG qua 
2 G 12<6 35 D 4:11 68 D 3-5 
3 C 6-3 | 36 Di U7 | 669 D 1-4 
4 C 10-9 | 37 oo | | ‘0 C 6-5 
5 C 9-10 38 C46 | 71 D 2-10 
6 Di 1-2 | 39 C 2-3 72 G lWs2 
7 ge 1-3 | 40 @ 12-7 73 C 7-10 
8 D 3-6 4] G 5s] | 74 D 24 
9 G T-7 42 G 12-2 | 75 Gi 2)? 

10 © 5=3 43 Co2e | 76 Dy 2-9 
11 me 1-16 44 % 77 C 5-9 
12 GC 2 45 S 78 CG 7-3 
13 @ Wf-3 46 ® 79 C 3-9 
14 D 3-8 47 Cc 80 C 9-6 
15 C wi-2 48 D 81 C 12-10 
16 Co9-5 CO 49 D 82 D 2-5 
17 C 1-8 50 C 83 G: 3.7 
18 C 10-7 51 C 84 C 1-5 
19 Ge a8 52 C 85 c 1-9 
20 Cle | 53 C 86 Cc 10-6 
21 C 8-1 54 D 87 G I<4 
22 G 1126 55 C 88 C 6.6 
23 i 343 56 D 89 Ci TOMI 
24 6 Gu 57 c 90 CG 7-5 
25 CG Ful 58 c 91 S 7x9 
26 D 4-9 59 D 92 C 3-5 
27 D 3«10 60 Cc 93 D 1-3 
28 D 4-10 61 C 94 Cc 9-8 
29 CG t=6 62 Cc 95 CG 4.4 
30 C 8.7 63 Cc 96 C 2-10 
31 G 1 64 c 97 C 902 
32 C 10-3 | 65 C 98 @ it<)1 
33 C c 99 G 1253 


4-2 | 66 


* Refers to specimen number in Test Plan. 
** Plate number - row - column. 
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TABLE 3 (Continued) 


Specimen Specimen | Specimen | 
Specimen Location Specimen Location Specimen Location 
No. * Code ** No. * Code ** | No. * Code ** 


100 G l6 126 3-4 | 152 4.5 
101 11-4 127 10.2 | 153 9-1 
102 4-10 128 5.10 | 154 4.5 
103 5-7 129 
104 9-9 130 
105 5-11 | 131 
106 9-11 132 
107 4-3 | 133 
108 10-4 | 134 
109 4-9 135 
110 12-8 | =: 136 
4-2 137 
207 138 
1-9 139 
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TABLE 81 


TEST RESULTS FOR B-58 DATA SET 
ORIGINAL STRAIN GAGE DATA 


2 
sequence A hasa@G value of 19.8 ksi’. 


2 


Specimen Sequence Failure Average No. 
Number Cycles Flights Cycles/ Flights 

13 A 564,937 1956 

14 A 782,822 2723 

15 A 767,827 2667 

16 A 559,334 1939 

17 A 1, 166,998 3710 

18 A 638,366 2221 
Average 746,714 2536 294 

19 C 1, 139,028 4047 

2] C 1,036,393 3679 

22 Cc 976,482 3466 

Re C 864,875 3048 

24 C 1, 153,512 4096 

25 C 1, 104, 352 3850 
Average 879,107 3697 237 

26 B 367,520 242] 

29 B 398,787 2598 

30 B 344,768 2260 

31 B 389,847 2554 

32 B 388,495 2542 

33 B 388,490 2542 a 
Average 379,651 2486 E52 

34 D 369,141 2511 

35 D 401,207 2717 

36 D 340,371 2338 

37 D 422,166 2887 

38 D 478,542 3211 

39 D 556,902 3786 —— 
Average 428,054 2908 147 


Sequence C same flight cycles as Spectrum A but ground (taxi) 


Sequence B 


cycles removed. 


flight refueling stress. 


hasagqy 


2 


value of 24.9 ksi’. 


Ground stress set equal to pre- 


Sequence D same flight cycles as Spectrum B but ground (taxi) 


cycles removed. 


flight refueling stress. 
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Ground stress set equal to pre- 


TABLE 82 


TEST RESULTS FOR F-106 DATA SET 
ORIGINAL STRAIN GAGE DATA 


Specimen Failure Average No. 
Number Sequence Cycles Flights Cycles/ Flight 
40 E 2,180,569 31,882 
4] E 2,855,495 41,743 
42 E 2,430,610 35,052 
43 E 2,850,984 41,703 
44 E 3,337,040 48,777 
45 E 2,161,772 31,595 
Average 2,636,078 38,458 68.5 
46 F 461,041 5,472 
47 F 465,350 5,530 
48 F 323,386 3,843 
49 F 704,095 8,371 
50 EF 987, 347 11, 781 
51 F 682, 153 8,092 
56 F 743,264 8,933 
57 F 631,791 7,480 _s 
Average 624, 803 7,441 83 


7 =: 6 
Sequence E has ag_ value of 18.23 x 10 
Sequence F has ag value of 41.90 x 10 


a is a weighted average which was calculated by the following equation 


N 2 
2 7 
2 ; 
o is i=l 
weighted N 
n. 
1 
i=l 
whe re 2 
O. = variance of flight i about the mean of the 
absolute stress peaks and valleys of flight 1 
n. = number of stress reversals in flight i 
N = number of flights in data set 
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Specimen 
Number Sequence 

75 A 

76 Folded 

vf See 

78 Paragraph 

79 Sikes 2 
Average 

107 A 

108 Bivariate 

109 Peak & Valley 

110 See 

lll Paragraph 

112 <a ee 
Average 


TABLE 83 


Failure 
Cycle 


414,810 
270,010 
320,866 
286,890 
302,445 
319,004 


471,162 
502,945 
540,842 
457,182 
497,683 
532,120 
500,322 
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Flights 


1712 
1127 
1335 
1197 
1262 
1326 
2084 
2219 
2417 
2018 
2123 
2325 
2198 


TEST RESULTS FOR SIMULATED B-58 SEQUENCE A 


Average No. 
Cycles/Flight 


240 


227 


TABLE 64 


TEST RESULTS FOR SIMULATED F-106 SEQUENCE E 


Specimen 
Number , sequence 
80 Hive visi 
81 Mean & 
82 Alternating 
83 See 
84 Paragraph 
85 5 20d 
Average 
86 Separate 
87 Peak & 
88 Valleys 
89 See 
91 Paragraph 
92 ee ees 
Average E 
93 Limited 
94 Separate 
95 Peak & 
96 Valley 
ot See 
98 Paragraph 
532.3 
Average 


Failure 
Flights 


Cycles 


536,792 
676,275 
595,466 
442,538 
683,415 


716,849 
608,555 


1,065,420 
1,410, 488 
799,311 
886,290 
661,348 
L, 1357, 803 
993,443 
485; 727 
1, 184,699 
923,274 
786,942 
837,690 
1,052,178 


928,418 
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18,931 
23,951 
21, jis 
14,616 
24, 187 
25,312 
21,352 


46, 169 
60,982 
34,667 
38,339 
26;:593 


49,204 
42,992 


35,505 
53,586 
41,752 
35,630 
37,852 
47,619 


42,002 


Average No. 
Cycles/Flight 


28.5 


A an 


22.1 


TABLE 85 
TEST RESULTS FOR SIMULATED B-58 SEQUENCE B, C, AND D 


133 


Specimen Failure Average No. 
Number Sequence Cycles Flights Cycles/ Flight 
125 C 556, 685 2487 
126 Bivariate 910, 042 4071 
127 Peak and 888, 343 3981 
128 Valley 620,371 2770 
130 See 581,368 2607 
13] Paragraphs 505,582 2268 
§.-l. 2 and 
5.4.1 —— 
Average 677, 065 3031 223 
132 B 264, 329 2161 
133 Bivariate 315,276 2547 
134 Peak and 360, 874 2959 
135 Valley 243,023 1998 
136 See 282,724 2295 
137 Paragraphs 269, 424 2201 
5.1.2 and 
5.4.1 . ss 
Average 289,275 2360 122 
140 D 414,844 3564 
14] Bivariate 499,617 4300 
142 Peak and 385, 938 3315 
143 Valley 397,527 3413 
See 
Paragraghs 
5.1.2 and 
5.401 | — 
Average 424,481 3648 116 


Specimen 
Number 


113 
114 
115 
116 
Lie 
118 


Sequence 


F 
Separate 
Peak & Valleys 
See 
Paragraphs 
5.2.2 and 
5.4.2 
Average 


TABLE 86 


TEST RESULTS FOR SIMULATED F-106 
SEQUENCE F 


Failure 
Cycles 


192,248 
209, 853 
308, 147 
175, 749 
158,982 
172,224 


202, 867 


134 


Flights 


6721 
7344 
10736 
6175 
S551 
6043 


7096 


Average No. 
Cycles/ Flight 


28.6 
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